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This  paper  introduces  important  new  (unctions  for  analytic  solution  ol  Lanch- 
ester-type  equations  of  modem  warfare  for  combat  between  two  homogeneous 
forces  modeled  by  power  attrition-rate  coefficients  with  ‘no  offset '  Tabulations 
of  these  Lanchester -Clifford-Schlafh  (or  LCS)  functions  allow  one  to  study  this 
particular  variable-coefficient  model  almost  as  easily  and  thoroughly  as  Lanch¬ 
ester  s  classic  constant-coefficient  one.  LCS  functions  allow  one  to  obtain 
important  information  (in  particular,  force-annihilation  prediction)  without  having 
to  spend  the  time  and  effort  of  computing  force-level  trajectories.  The  choice  ol 
these  particular  functions  is  based  on  theoretical  considerations  that  apply  in 
general  to  Lanchester-type  equations  of  modem  warfare  and  provide  guidance 
for  developing  other  canonical  functions.  Moreover,  our  new  LCS  functions  also 
provide  valuable  information  about  related  variable-coefficient  models  Also,  we 
introduce  an  important  transformation  of  the  battle  s  time  scale  that  not  only 
simplifies  the  force-level  equations,  but  also  shows  that  relative  fire  effectiveness 
and  intensity  of  combat  are  the  only  two  weapon-system  parameters  determin¬ 
ing  the  course  of  such  variable-coefficient  Lanchester-type  combat 


IN  AN  EARLIER  paper  (Taylor  and  Brown  [1976]),  we  showed  how 
to  solve  variable-coefficient  Lanchester-type  equations  of  modern 
warfare  for  combat  between  two  homogeneous  forces.  In  that  paper,  we 
introduced  canonical  hyperbolic-like  Lanchester  functions  for  construct¬ 
ing  the  solution.  Unfortunately,  with  only  these  previous  results  one  is 
limited  to  computing  force-level  trajectories  and  cannot  gain  a  real 
understanding  of  qualitative  model  behavior  (e.g.  force  annihilation) 
without  extensive  numerical  computations  (and  only  then  for  specific 
values  of  model  parameters).  Since  the  appearance  of  our  earlier  work, 
several  mathematical  discoveries  (Taylor  and  Comstock  [1977],  Taylor 
[1979b])  have  provided  new  qualitative  insight  about  the  behavior  of  this 
combat  model.  We  wish  to  show  here  how  these  new  results  allow 
parametric  analysis  of  combat  modeled  by  power  attrition-rate  coeffi¬ 
cients  with  somewhat  the  same  facility  as  allowed  by  F.  W.  Lanchester’s 
classic  constant-coefficient  model.  In  order  to  obtain  this  analysis  capa- 
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bility,  however,  one  must  redefine  the  Lanchester-Clifford-Schlafli  (or 
LCS)  functions,  which  we  introduced  in  Taylor  and  Brown  [1976]. 

It  is  important  for  the  military  OR  analyst  to  have  a  clear  understand¬ 
ing  of  how  the  initial  force  ratio  and  weapon-system-capability  parame¬ 
ters  interact  to  determine  a  battle’s  outcome.  One  is  consequently  inter¬ 
ested  in  developing  insights  into  the  dynamics  of  combat  by  explicitly 
portraying  the  relation  between  the  various  factors  of  the  combat- 
attrition  process  and  battle  outcome.  Modeling  battle  termination  is  a 
somewhat  controversial  topic  (Taylor  [1979a]),  and  no  mathematical 
theory  exists  for  other  than  determining  zero  points  of  solutions 
(i.e.  force  annihilation)  to  such  differential-equation  models  (Taylor 
[1979b]).  However,  it  is  of  considerable  utility  just  to  be  able  to  easily 
predict  the  occurrence  of  force  annihilation  in  simulated  Lanchester- 
type  combat.  One  is  always  interested  in  determining  what  conditions 
lead  to  the  annihilation  of  an  enemy  force,  since  such  an  occurrence  (of 
course)  guarantees  victory.  Although  actual  battles  rarely  go  completely 
to  annihilation,  a  commander  may  decide  to  terminate  an  engagement 
once  he  anticipates  that  annihilation  is  possible,  and  hence  force-anni¬ 
hilation  conditions  may  be  useful  in  modeling  engagement  termination. 
Additionally,  a  commander  would  seek  to  avoid  engagements  in  which 
his  own  force  could  be  annihilated,  and  such  conditions  may  provide,  for 
example,  valuable  information  for  the  modeling  of  engagement  avoidance. 

In  our  earlier  paper  (Taylor  and  Brown  [1976]),  we  gave  various 
examples  of  hyperbolic-like  Lanchester  functions  (in  particular,  the  LCS 
functions,  which  arise  from  power  attrition -rate  coefficients  with  “no 
offset”).  Subsequent  research  by  Taylor  and  Comstock  has  revealed, 
however,  that  these  canonical  LCS  functions  must  be  redefined  to  permit 
force-annihilation  prediction  from  initial  conditions  without  having  to 
spend  the  time  and  effort  to  compute  force-level  trajectories.  It  then 
became  obvious  that  the  entire  topic  of  representing  the  solution  to  such 
Lanchester-type  equations  in  terms  of  general  Lanchester  functions 
(GLF )  should  be  critically  reexamined.  Consequently,  we  developed  new 
general  considerations  for  the  selection  of  canonical  Lanchester  functions 
(Taylor  and  Brown  [1977a]).  Based  on  these  considerations,  we  also 
developed  new  LCS  functions  for  the  special  case  of  power  attrition -rate 
coefficients  with  “no  offset”  (modeling,  for  example,  weapon  systems 
with  the  same  maximum  effective  range)  which  are  presented  here.  These 
power  Lanchester  (i.e.  LCS)  functions  are  significant,  not  only  because 
they  correspond  to  attrition-rate  coefficients  modeling  a  large  class  of 
combat  situations  of  interest,  but  also  because  they  yield  valuable  infor¬ 
mation  about  other  related  canonical  Lanchester  functions,  e.g.  the  offset 
power  Lanchester  functions  (see  Taylor  and  Brown  [1978]  and  Section 
10  below).  With  the  availability  of  tabulations  of  these  new  LCS  func- 
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tions,  one  can  study  this  model  almost  as  easily  and  thoroughly  as 
Lanchester’s  classic  constant -coefficient  one.  Such  models  are  important 
for  developing  insights  into  the  dynamics  of  combat  (Bonder  and  Honig 
[1971],  Taylor  [1980a]). 

The  results  of  this  paper  are  also  important  for  understanding  complex 
operational  differential-equation  models  that  are  widely  used  in  both  the 
United  States  and  also  NATO  countries  as  defense-planning  tools  (see 
Huber  et  al.  [1975, 1979],  Taylor  [1979a]).  The  modern  high-speed,  large- 
scale  digital  computer  has  made  it  possible  to  develop  and  use  such 
complex  Lanchester-type  combat  models  (e.g.  see  Bonder  and  Farrell 
[1970];  Bonder  and  Honig;  Command  and  Control  Technical  Center 
(CCTC)  [1979]).  Nevertheless,  a  simple  combat  model  such  as  we  con¬ 
sider  here  may  yield  a  clearer  understanding  of  important  relations  that 
are  difficult  to  perceive  in  a  more  complex  model,  and  such  insights  can 
provide  valuable  guidance  for  subsequent  higher-resolution  computerized 
investigations.  As  Geoffrion  [1976]  has  emphasized,  one  can  use  a  sim¬ 
plified  auxiliary  model  for  understanding  the  basic  dynamics  and  behavior 
of  a  large-scale  complex  operational  model.  Furthermore,  one  can  fit  an 
analytical  model  to  data  generated  from  a  detailed  combat  simulation, 
and  thus  a  simple  analytical  model  like  the  one  considered  here  may 
provide  an  economical  framework  for  summarizing  simulation  output 
data  (see  Ignall  et  al.  [1978]  for  a  lucid  discussion  of  this  modeling 
strategy  in  a  nonmilitary  context). 

1.  VARIABLE-COEFFICIENT  LANCHESTER-TYPE  EQUATIONS  OF 
MODERN  WARFARE 

We  consider  the  following  variable -coefficient  Lanchester-type  equa¬ 
tions  of  modern  warfare  for  combat  between  two  homogeneous  forces  for 
x  and  y  >  0  (see  p.  45  of  Taylor  and  Brown  [  1976]  for  further  discussion) 


dx/dt  =  -a(t)y  with  x(0)  =  Xo, 

dy/dt  =  —b(t)x  with  y(0)  =  y0 , 


(1.1) 


where  t  =  0  denotes  the  time  at  which  the  battle  begins,  x(£)  and  yft) 
denote  the  numbers  of  X  and  Y  at  time  t,  and  aft)  and  bit)  denote  time- 
dependent  Lanchester  attrition-rate  coefficients,  which  do  not  explicitly 
depend  on  x  and  y.  In  particular,  both  a  and  b  depend  explicitly  upon 
time  (perhaps  via  an  intermediate  variable  such  as  range  r(f )),  but  a  does 
not  directly  depend  on  the  number  of  targets  x.  Although  combat  between 
two  military  forces  is  a  complex  random  process,  such  a  deterministic 
model  of  the  combat  attrition  process  is  frequently  employed  to  provide 
insights  into  the  dynamics  of  combat  (e.g.  see  Weiss  [1957];  Bonder  and 
Farrell;  Bonder  and  Honig;  Taylor  and  Parry  [1975],  Taylor  [1980a]). 
Moreover,  current  large-scale  operational  models  (e.g.  see  Bonder  and 
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Farrell,  Bonder  and  Honig,  CCTC)  more  or  less  take  (1.1)  as  the  point 
of  departure  for  their  development  through  the  process  of  model  enrich¬ 
ment  (see  Morris  [1967]  for  a  lucid  discussion  of  this  enrichment  process). 
For  example,  in  the  detailed  VECTOR-2  operational  model  (e.g.  see 
CCTC),  the  attrition -rate  coefficients  are  nonautonomous  and  depend 
(in  quite  a  complicated  fashion)  on,  not  only  the  engagement  conditions 
(e.g.  range  between  firer  and  target,  target  and/or  firer  motion,  posture, 
etc.),  but  also  the  number  of  firers  and  targets. 

Equations  1.1  are  usually  taken  to  model  combat  in  which  both  sides 
use  aimed  fire  and  target  acquisition  times  are  independent  of  the  number 
of  firers  and  targets  (see  Taylor  [1974,  1980a],  Taylor  and  Brown  [1976] 
for  further  details).  Other  forms  of  Lanchester-type  equations  appear  in 
the  literature,  but  we  will  not  consider  them  here  (see  Dolansky  [1964], 
Taylor  [1974,  1979a,  1980a]).  The  Lanchester  attrition-rate  coefficients 
a(t )  and  b(t)  depend  on  such  variables  as  force  separation,  tactical  posture 
of  targets,  rate  of  target  acquisition,  firing  doctrine,  firing  rate,  and  so 
forth  (e.g.  see  Bonder  [1965,  1967,  1970];  Bonder  and  Farrell).  Bonder 
[1965]  (see  also  Bonder  and  Farrell)  has  stressed  the  importance  for 
evaluating  weapon  systems  of  such  variable  coefficient  differential  com¬ 
bat  models  to  represent  temporal  variations  in  firepower  on  the  battle¬ 
field. 

We  assume  that  a(t)  and  b(t)  are  defined,  positive,  and  continuous  for 
t0  <  t  <  +»  with  to  <  0  (see  Taylor  and  Brown  [1976],  Taylor  [1979b] 
for  further  discussion).  We  further  assume  that  a(t)  and  6(f)  are  such 
that  their  right-hand  limits  exist  at  f0,  with  +oo  allowed  as  a  possibility: 
we  define  a(f0)  as  limj^+aff)  and  similarly  define  6(fo).  Note  that  the 
values  0  and  +oo  are  possible  for  a(t)  and  6(f)  only  at  t  =  to.  For 
convenience,  we  introduce  the  notation  a(t)  E  L(to,  T)  to  mean  that 
Jj  a(t)dt  exists.  We  also  assume  that  a(f)  and  6(f)  E  L(f0,  T)  for  any 
finite  T  2r  fo.  It  follows  that,  for  example,  a(t)  E  L(to,  +°°)  implies  that 
limr— +«  /to  a(t)dt  =  +oo.  We  will  further  take  a(t)  and  6(f)  to  be  given 
in  the  form  a(t)  =  Kg(t)  and  6(f)  =  6*6(f),  where  K  and  kb  are  positive 
constants  chosen  so  that  a(f)/6(f)  =  kjko  if  and  only  if  g(t)  =  6(f)  (see 
Taylor  [1979b,  1982]).  In  other  words,  ka  and  kb  are  basically  “scale 
factors,”  which  are  useful  for  parametric  study  of  battle  outcomes  as 
related  to  various  system  parameters.  This  factorization  of  a(t)  and  6(f) 
is  not  used  directly  in  (1.1),  but  is  implicit  in  constructing  the  general 
Lanchester  functions  used  to  represent  the  analytical  solution  to  (1.1) 
(see  pp.  441  and  448  of  Taylor  [1979b]).  It  is  also  convenient  to  introduce 
the  combat -intensity  parameter  and  the  relative -fire-effectiveness  pa¬ 
rameter  Xfi  defined  by 

=  y/kokt,,  and  \r  =  kjkt,.  (1.2) 

(See  Taylor  and  Brown  [1978],  Taylor  [1979a,  1980b)  for  further  details.) 
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The  X  force  level  as  a  function  of  time,  x(t),  may  be  represented  as 
(Taylor  and  Brown  [1976]) 

x(t)  =  Xo|Cy(0)C*<t)  -  smsdt) |  (L8) 

-  yo'/xi?|Cx(0)Sx(f)  —  Sx(0)Cx(O). 

where  the  hyperbolic -like  general  Lanchester  functions  (GLF)  Cxit)  and 
Sx{t )  are  linearly  independent  solutions  to  the  X  force-level  equation 

d2x/dt 2  -  \[\/a(t)]da/dt\dx/dt  -  a(t)b(t)x  =  0,  (1.4) 

with  initial  conditions  Cx(to )  =  1,  \l/a(to)\dCx/dt{to)  =  0,  SxUo)  =  0,  and 
|l/a(fo)|dSx/df(to)  =  1  When  (for  example)  a(fo)  =  0  or  +a>,  an 
initial  value  such  as  |l/a(to)|dCx/df(f0)  should  be  interpreted  as 
lim,_^+[|l/aU)|dCx/dl(t)].  Taylor  and  Comstock,  and  Taylor  [1979b] 
introduced  and  studied  exponential -like  GLF.  In  Taylor  and  Brown 
[1976],  we  have  discussed  the  representation  of  force  levels  in  terms  of 
GLF  and  have  shown  that  these  two  types  (i.e.  hyperbolic- like  and 
exponential -like  GLF)  are  essentially  the  only  kinds  of  GLF,  but  that 
the  hyperbolic-like  ones  are  to  be  preferred. 

2.  GENERAL  FORCE-ANNIHILATION-PREDICTION  CONDITIONS 

The  following  theorem  generalizes  Lanchester’s  famous  square  law  to 
variable-coefficient  combat  (see  Taylor  [1979b]  for  proof  of  a  more 
general  result). 

THEOREM  1  (Taylor  and  Comstock).  Assume  that  either  a(t)  £  L(0,  +«) 
or  b(t)  (£  L(0,  +oo).  Then  the  X  force  will  be  annihilated  in  finite  time  if 
and  only  if 

xo/yo  <  ^|[Cx(0)  -  Q*Sx(0)]/[Q*Cy(0)  -  Sy(0)]),  (2.1) 

where  the  parity -condition  parameter  Q*  is  unique  and  given  by 

lim^.[Sx(t)/Cx(t)]  =  1  fOT  =  l/llim^SyttVCytt)]).  (2.2) 

An  answer  to  the  seemingly  simple  question  “Who  will  be  annihilated 
in  battle?”  requires  a  significant  extension  of  the  theory  of  the  real  zeros 
of  nonoscillatory  (in  the  strict  sense)  solutions  to  the  general  second- 
order  linear  differential  equation  (Taylor  [1979b]).  Furthermore,  consid¬ 
eration  of  Theorem  1  shows  that  the  power  Lanchester  (or  LCS)  func¬ 
tions  introduced  in  Taylor  and  Brown  [1976]  were  inappropriately  de¬ 
fined  (see  Taylor  and  Brown  [1977a]  for  further  details).  It  is,  therefore, 
the  purpose  of  this  paper  to  appropriately  redefine  the  power  Lanchester 
functions  in  light  of  Theorem  1. 
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3.  COMBAT  MODELED  WITH  POWER  ATTRITION-RATE 
COEFFICIENTS 

A  large  class  of  tactical  situations  of  interest  can  be  modeled  with  the 
following  general  power  attrition-rate  coefficients 

a(t)  =  K,(t  +  KsY,  and  b(t)  =  +  Ks  +  K0Y,  (3.1) 

where  fi,  v,  Ks,  Kn  ^  0.  Taylor  and  Brown  [1976]  discuss  the  modeling 
roles  of  K$  and  Ko .  We  will  call  Ks  the  starting  parameter,  since  it  allows 
us  to  model  battles  that  begin  within  the  maximum  effective  ranges  of 
both  opponents.  We  will  call  Kn  the  offset  parameter,  since  it  allows  us 
to  model  battles  between  opposing  weapon  systems  with  different  maxi¬ 
mum  effective  ranges.  We  also  observe  that  t0  =  —Ks- 

The  above  nomenclature  is  motivated  by  Bonder’s  [1965]  model  of  a 
constant -speed  attack  against  a  static  defensive  position 

dx/dt  =  —ct(r)y,  and  dy/dt  =  —0(r)x,  (3.2) 

where  r  denotes  the  range  between  opposing  forces,  and  «(r)  and  0(r) 
denote  range-dependent  attrition -rate  coefficients.  Range  is  related  to 
time  by 

r(t)  =  r0  —  vt,  (3.3) 

where  0  <  t  <  r0/v,  r0  denotes  the  opening  range  of  battle,  and  v  >  0 
denotes  the  constant  attack  speed.  For  example,  consider  the  constant- 
speed  attack  of  a  homogeneous  Y  force  against  the  static  defensive 
position  of  a  homogeneous  X  force  (see  Figure  1).  The  basic  idea  is  that 
force  separation  (i.e.  range  between  the  opposing  forces)  changes  over 
time  and  that  the  fire  effectiveness  of  (for  example)  a  single  Y  firer, 
denoted  as  a(r),  depends  on  the  force  separation. 

In  many  cases  of  tactical  interest,  we  may  model  the  fire  effectiveness 
of  the  y  weapon  system  as  a  function  of  range  with  (see  pp.  196-200  of 
Bonder  and  Farrell) 

[ao(l  -  r/rJY  for  0  <  r  <  ra, 
at(r)  =  \  (3.4) 

l  0  for  r„  <  r, 

where  r„  denotes  the  maximum  effective  range  of  the  Y  weapon  system 
and  n  >  0  models  the  range  dependency  of  Y' s  attrition -rate  coefficient 
(see  Figure  2).  We  model  0(r)  similarly,  with  corresponding  quantities  rB 
and  v. 

Substituting  (3.3)  and  (3.4)  into  (3.2),  we  find  that  K0  =  ( rB  -  ra)/u 
and  Ks  =  (r„  -  r0)/v,  and  that  ka  =  ao(u/rY  and  kb  =  Mv/rg)".  Thus,  Ko, 
Ks  >  0  if  and  only  if  rB  >  r„  >  r0.  Moreover,  for  this  particular  application 
(and  this  situation  is  typical),  the  attrition-rate  coefficients  are  techni- 
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cally  not  defined  in  the  tactical  scenario  for  t  >  tmtt.  Nevertheless,  one 
can  conceptually  embed  Bonder’s  tactical  model  with  these  time-depend¬ 
ent  attrition-rate  coefficients  in  the  mathematical  model  in  which  (3.1) 
are  assumed  to  hold  for  all  t  >  0  as  long  as  one  is  careful  not  to  use  any 
quantities  computed  from  the  mathematical  model  outside  range  of 
definition  of  the  tactical  scenario.  In  particular,  one  must  verify  that 
force  annihilation  occurs  within  the  range  of  definition  of  the  tactical 
scenario  (e.g.  before  t, m  for  Bonder’s  tactical  model),  and  this  require¬ 
ment  generates  the  need  of  computing  the  time  at  which  force  annihila¬ 
tion  occurs  in  the  mathematical  model  in  which  the  coefficients  are 
assumed  to  hold  for  all  t  >  0.  We  raise  this  point  again  in  Sections  6  and 
8  and  illustrate  it  in  the  examples  given  in  Section  9. 


Figure  1.  Diagram  of  Bonder’s  constant-speed  attack  model.  Force 
separation,  r(t),  is  given  by  r(t)  =  r0  -  vt. 


When  the  offset  parameter  is  equal  to  zero  (i.e.  K0  =  0),  the  coefficients 
(3.1)  reduce  to  the  following  power  attrition- rate  coefficients  with  “no 
offset  “ 

a(t)  =  +  Ks)“,  and  b(t)  =  kbU  +  Ks)\  (3.5) 

As  we  have  just  seen  above  in  Bonder’s  model,  these  coefficients  model 
(for  example)  combat  between  weapon  systems  with  the  same  maximum 
effective  range,  so  that  there  is  no  “offset”  in  the  capabilities  of  the 
opposing  systems  to  “reach  out”  on  the  battlefield.  It  is  the  purpose  of 
this  paper  to  introduce  new  power  Lanchester  functions  that  facilitate 
force-annihilation  prediction  (and  also  determination  of  how  long  the 
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battle  will  last)  for  “aimed-fire”  combat  modeled  by  the  power  attrition- 
rate  coefficients  with  “no  offset”  (3.5).  The  results  of  this  section  show 
how  the  physical  characteristics  of  the  weapon  systems  and  environment 
are  related  to  these  coefficients. 

4.  A  TRANSFORMATION  TO  NORMALIZE  THE  BATTLE’S  TIME 
SCALE  BY  THE  INTENSITY  OF  COMBAT 

In  this  section  we  show  how  transformation  of  the  battle’s  time  scale 


Figure  2.  Dependence  of  Y's  attrition-rate  coefficient  a(r)  on  the 
exponent  n  with  the  maximum  effective  range  of  the  weapon  system  and 
kill  rate  at  zero  range  held  constant.  [Notes:  (1)  The  maximum  effective 
range  of  the  system  is  denoted  as  r„  =  2000  meters.  (2)  a(0)  =  t»o  =  0.6X 
casualties/(unit  time  X  number  of  Y  firers)  denotes  the  weapon-system 
W ill  rate  for  Y  at  zero  force  separation  (range).  (3)  The  opening  range  of 
battle  is  denoted  as  r0  =  1250  meters  and  (as  shown)  r0  <  r„.) 

provides  important  insight  into  the  parametric  dependence  of  the  course 
of  combat.  Accordingly,  we  introduce  the  new  independent  variable  t 
defined  by  (see  Taylor  and  Brown  [1976,  1977a],  Taylor  [1979b]  for 
further  details) 

t(0  =  f  < Ja(s)b(s)ds ,  (4.1) 

Jto 

and  let  r0  denote  r(0).  As  is  readily  seen,  this  transformation  is  well 
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defined  and  invertible.  We  observe  that  to  <  0  implies  that  r0  >  0.  If  we 
denote  the  “average  intensity  of  combat”  as  Va(t)6(t),  then 


=  |  (1/t)  I  Va(s)6(s)ds  ji  =  r  —  t0.  (4.2) 


The  substitution  (4.1)  transforms  (1.4)  into 


d2x/dr2  -  m\d  In  R(t)/dr\dx/dr  -  x  =  0,  (4.3) 

with  initial  conditions  x(r0)  =  Xo  and  |l/Vfl(0) \dx/dT(r0)  =  -yo,  where 
R(t)  =  a(t)/b(t).  Equation  4.3  is  highly  significant  (see  Section  10  below) 
because  it  clearly  shows  that  the  course  of  combat  depends  on  just  the 
two  weapon-system  parameters:  (1)  R(t)  =  a{t)/b(t),  the  relative  fire 
effectiveness  (V  to  X)  of  the  opposing  weapon-system  types,  and  (2)  I(t) 
=  y/a(t)b(t),  the  intensity  of  combat  (through  (4.1),  which  relates  I(t)  to 
r).  In  particular,  from  (4.3)  we  see  that  the  nature  of  temporal  variations 
in  relative  fire  effectiveness  will  have  a  significant  effect  on  the  course 
of  combat  (see  Taylor  [1980b]  for  further  details). 

For  the  power  attrition-rate  coefficients  with  no  offset  (3.5),  the 
transformed  X  force-level  equation  becomes 

d2x/dr 2  +  |(2 q  -  l)/r|dx/dT  -  x  =  0,  (4.4) 

with  initial  conditions  x(r0)  =  *o  and  \(T/2)*-ldx/dT\„'  =  -y0>/A*|A//(/i 
+  v  +  2)|.  Here 

q  =  («/  +  l)/(n  +  v+2),  (4.5) 

and 


r(t)  =  |2A7/(ji  +  ^  +  2)|(t  +  Ks)(^2>/2.  (4.6) 

Let  us  observe  that  0  <  q  <  1  when  n  and  v  >  —  1.  Furthermore,  q  >  lA  if 
and  only  if  dR/dt  <  0,  i.e.  R(t)  is  a  strictly  decreasing  function  of  time. 


5.  LANCHESTER-CLIFFORD-SCHLAFLI  (LCS)  FUNCTIONS 

Consider  the  function  Fa(x)  defined  by  the  power  series 

Fa(x)  =  r(a)  EL.  (x/2)»/\k\T(k  +  a)|.  (5.1) 

For  a  ¥“  0,  —1,  —2,  —  ,  the  radius  of  convergence  for  F„(x)  is  infinite  by 
the  ratio  test  for  convergence  of  power  series  (Knopp  [1956]).  Hence, 
Fa(z)  is  an  entire  function  of  the  complex  variable  z  =  x  +  iy  with  an 
essential  singularity  at  the  point  of  infinity.  Now  consider  the  function 
HJx )  defined  by  the  infinite  series 

hjx)  =  r(«)  Eu  (x/2)«*-yi*!r(*  +  «  +  dj.  (5.2) 

Observing  that  Ha(x)  =  (l/a)(x/2)2"ForM(x),  we  see  that  for  a  >  0  the 
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infinite  series  (5.2)  is  uniformly  convergent  on  compact  subsets  of  the 
complex  plane.  One  can  also  readily  deduce  the  recursive  relation  F„(x) 
=  Fo+i(x)  +  t(x/2)2/[a(a  +  l)|F0+2(x).  We  will  call  the  functions  F„(x) 
and  H„(x)  Lanchester-Clifford-Schlofli  (LCS)  functions.  Other  properties 
are  readily  deduced  and  are  given  in  Table  I. 

Although  the  solution  of  the  X  force-level  Equation  1.4  with  the  power 
attrition-rate  coefficients  (3.5)  may  be  expressed  in  terms  of  known 
higher  transcendental  functions  (see  Taylor  [1974],  Taylor  and  Brown 
[1976],  Taylor  and  Comstock),  we  have  chosen  to  introduce  the  LCS 
functions,  since  tabulations  of  these  other  functions  are  only  available 
for  a  very  restrictive  range  of  parameter  values  of  interest  in  Lanchester 
combat  theory.  For  example,  we  can  construct  such  solutions  with 
modified  Bessel  functions  of  the  first  kind  of  fractional  order,  but 
tabulations  of  these  (e.g.  see  Abramowitz  and  Stegun  [1964])  exist  only 
for  a  restrictive  set  of  values  of  the  order  p  (i.e.  p  =  ±‘/«,  +V»,  ±'A,  ±%,  ± 


TABLE  I 

Properties  op  the  LCS  Functions  F.{x)  and  H.(x) 

1.  dFJdx  -  (i/2)1J*W.(i). 

2.  dHJdx  =  (x/2)^'F.(x). 

3.  F.(0)  -  1. 

4.  HJ 0)  -  0  for  a  >  0. 

5.  F.(x)F,..{x )  —  H.(x)Hi^,(x)  =  1  for  all  x  where  a  is  neither  an  integer  nor  zero. 

6.  dF./dx(0)  -  0. 

7.  l(x/2),udHJdx(x)l^,=  l. 

8.  F,n(x)  =  coah  x. 

9.  Hifiix)  =  sinh  x. 


3/4 ),  where  p  =  (m  +  1)/(m  +  v  +  2).  Furthermore,  there  are  no  tabulations 
of  functions  corresponding  to  the  quotient  of  two  GLF.  Consequently, 
we  have  introduced  our  new  LCS  functions,  which  provide  much  of  the 
information  desired  about  such  battles.  The  naming  of  our  LCS  functions 
follows  from  the  facts  that  a  function  similar  to  F„(x)  was  introduced  by 
Ludwig  Schlafli  [1867/68]  and  that  a  related  one  appears  in  a  posthumous 
fragment  of  the  great  English  geometer  William  Kingdon  Clifford 
[1882]. 

The  function  F„(x)  satisfies  the  linear  second-order  ordinary  differen¬ 
tial  equation 

d2FJdx2  +  |(2 a  -  l)/x|dF0/dx  -  Fa  =  0,  (5.3) 

with  initial  conditions  F„(0)  =  1  and  dFJdx( 0)  =  0,  while  Ha(x)  satisfies 
d2HJdx2  -  | (2«  -  l)/x] dHJdx  -  Ha  =  0,  (5.4) 

with  initial  conditions  (for  a  >  0)  Ha( 0)  =  0  and  |(x/2),-2"d//„/dx(x))I_o 
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=  1.  Thus,  |Fa,  H i_a)  is  a  fundamental  system  of  solutions  to 

d2F/dx2  +  { (2a  -  1  )/x\dF/dx  -  F  =  0,  (5.5) 

with  Wronskian  W(Fa,  Hx-a)  =  (x/2)1-2".  It  follows  that  the  GLF  for  the 
X  and  V  force-level  equations  for  combat  modeled  by  (1.1)  with  attrition- 
rate  coefficients  (3.5)  are  given  by 

CxU )  =  Fv(r),  Sx(t)  =  |X,/(M  +  *  +  2)\2*~lHp(r),  (5.6) 

CAt)  =  Fp(r),  SAt)  =  |X,/(M  +  *  +  2)|,-*'//,(t),  (5.7) 

where  q  is  given  by  (4.5),  p  =  1  -  q,  and  r(t)  is  given  by  (4.6). 

If  we  define 

T„(x)  =  Hl-a(x)/Fa(x)  (5.8) 

and  let  Tx(t)  =  Sx(t)/CAt)  denote  a  hyperbolic-like  GLF  corresponding 
to  the  hyperbolic  tangent,  then 

TAt)  =  \\,/(n  +  *  +  2)}*-'r„(r).  (5.9) 

Taylor  (1979b)  shows  that  for  n  and  v  >  —1,  Tx(t)  is  a  strictly  increasing 
function  with  range  [0,  1/Q*]  for  t  E  (0,  -H»)  and  that 

Q*  =  |r(q)/r(p)||X//(*i  +  ,  +  2)|1-2^.  (5.io) 

Consequently  (2.2)  and  (5.9)  yield  that  T„(x)  is  strictly  increasing  with 
T„(0)  =  0  and 

lim^TAx)  =  T(1  -  a)/r(a).  (5.11) 

6.  USE  OF  LCS  FUNCTIONS  FOR  ANALYZING  COMBAT 

The  Lanchester-Clifford-Schlafli  (LCS)  functions  Fa(x)  and  H„(x)  are 
very  useful  for  analyzing  “aimed-fire”  combat  modeled  by  the  power 
attrition -rate  coefficients  with  “no  offset”  (3.5).  Here,  we  assume  that 
the  attrition -rate  coefficients  (3.5)  hold  for  all  t  >  0  in  the  mathematical 
model  (1.1).  Recall  (see  Section  3  above)  that  one  must  be  careful  not  to 
use  any  results  computed  from  the  mathematical  model  out  of  the  range 
of  definition  of  the  tactical  scenario  describing  the  tactical  situation 
considered  in  any  particular  application.  For  such  combat,  these  LCS 
functions  may  be  used  to  (1)  compute  force  levels  as  a  function  of  time, 
(2)  predict  force  annihilation,  and  (3)  compute  the  time  of  force  annihi¬ 
lation.  We  will  now  show  how  to  obtain  this  information. 

According  to  (1.3),  (5.6),  and  (5.7),  we  may  write  the  X  force  level  as 

x(t)  =  Xo|Fp(ro)F,(T)  -  Hv(to)Hp(t)\ 

-  yo'/\K\\i/(n  +  »  +  2))2*-1|F,(7-o)//,(t)  -  Hp(ro)Fq(r)). 


(6.1) 
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From  (5.11)  and  (6.1)  (see  Taylor  [1979b]  for  details),  we  may  conclude 
the  following  force-annihilation-prediction  result.  (Alternatively,  we  may 
substitute  (5.6),  (5.7),  and  (5.10)  into  Theorem  1  to  obtain  Theorem  2.) 

THEOREM  2:  Consider  combat  between  two  homogeneous  forces  modeled 
for  all  t  >  0  by  (1.1)  with  power  attrition-rate  coefficients  (3.5).  Assume 
that  n  and  v  >  -1.  Then  the  X  force  will  be  annihilated  in  finite  time  if 
and  only  if 

*o/yo  <  VMX,/(/«  +  r  +  2)r',(r(p)/r(g))|[F,(T0)  (6 

-  [r(9)/r(p)]Hp(ro)]/lFp(ro)  -  [r(p)/r(9)]H,(r0)]|. 

When  T0  —  0  (i.e.  Ks  -  0),  the  X  force  will  be  annihilated  to  finite  time  if 
and  only  if 

xo/yo  <  VM A,/(m  +  v  +  2))9~p(r(p)/T(q)).  (6.2a) 

However,  one  must  verify  that  force  annihilation  does  not  occur  out  of 
the  range  of  definition  of  the  tactical  scenario  for  any  particular  appli¬ 
cation  (e.g.  after  t for  Bonder’s  constant-speed  attack  model  considered 
in  Section  3  above).  Turning  now  to  the  determination  of  the  time  at 
which  annihilation  occurs  for  the  mathematical  model  in  which  the 
attrition -rate  coefficients  (3.5)  have  been  assumed  to  hold  for  all  t  >  0, 
we  see  that  when  (6.2)  is  satisfied,  the  time  to  annihilate  the  X  force,  tax, 
is  determined  by  x(tax)  =  0.  It  follows  that 

Tv[r(tx)]  =  [x<,Fp(to)  +  yo>/MX,/(M  +  »  +  2)FP  (g  3) 

■Hp(To)]/[x0Hq(To)  +  yoVx^[Xr/(M  +  »  +  2)rpF,(ro)l, 
or,  more  explicitly, 

tx  -  r-,|77l[[x0F,P(ro)  +  y0V\H[\,/(n  +  r  +  2)]"  (g 

•Hp(ro)]/lxoH,(ro)  +  y0Vx^[X,/(M  +  r  +  2)rPF,(r0)]]}, 
where  r-*  and  T^x  denote  inverse  functions. 

7.  AVAILABLE  TABULATIONS  OF  LCS  FUNCTIONS 

Tabulations  of  the  Lanchester-Clifford-Schlafli  functions,  which  are 
given  in  Taylor  and  Brown  [1977b,  1977c],  are  available  from  the  Na¬ 
tional  Technical  Information  Service.  These  reports  contain  five-deci- 
mal-place  tables  of  the  hyperbolic-like  LCS  functions  F„(x),  and 

T„(x)  for  values  of  the  argument  x  =  0.00(0.01)2.00(0.1)10.0  and  various 
values  of  the  order  a.  The  short  table  (Taylor  and  Brown  [1977c]) 
contains  tabulations  for  eleven  values  of  a  in  the  range  (0,  1)  correspond¬ 
ing  to  m,  *  —  0,  1,  2,  3;  while  the  longer  table  (Taylor  and  Brown 
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[1977b])  contains  tabulations  for  26  values  of  a  in  the  same  range 
corresponding  to  n,  v  =  0,  V*,  'A,  1,  llA,  2,  3.  As  we  have  seen  above  in 
Section  3  (see  (3.2),  (3.4),  and  Figure  2),  such  values  of  n  and  v  allow  one 
to  analyze,  for  example,  a  wide  variety  of  range  capabilities  for  weapon 
systems  in  Bonder’s  [1965]  constant -speed  attack  model  (3.2).  These 
tables  have  been  calculated  by  the  recursive  methods  given  in  Section  8 
of  Taylor  and  Brown  [1976]. 

8.  OUTLINE  OF  COMPUTATIONAL  PROCEDURE 

The  above-mentioned  tabulations  of  these  new  LCS  functions  make 
the  analysis  of  an  important  class  of  Lanchester-type  battles  a  compar¬ 
atively  easy  matter.  Before  we  consider  numerical  examples  to  show  that 
insights  may  be  easily  obtained  into  the  dynamics  of  combat,  let  us 
outline  the  general  computational  procedure  (based  on  the  results  given 
in  Section  6)  that  one  follows  in  the  analysis  of  such  combat.  Accordingly, 
the  basic  steps  involved  are  as  follows: 

(1)  Determine  from  (6.2)  whether  the  X  force  can  be  annihilated, 

(2)  If  annihilation  is  possible,  determine  the  time  of  the  X  force’s 

annihilation  as  follows: 

(a)  Compute  Tq(rx)  by  (6.3)  [here  rax  =  r(t0*)], 

(b)  Using  interpolation,  determine  rax  from  the  appropriate  tab¬ 
ulation  of  Tq,  and 

(c)  Using  (4.5),  compute  tx  =  r-,( rx). 

Note  from  the  above  that  these  two  determinations  involve  only  the 
initial  force  ratio  Uo  =  Xo/yo  (and  not  the  individual  initial  force  levels 
themselves).  Additionally,  one  must  verify  that  such  numerical  results 
hold  within  the  range  of  definition  of  the  tactical  scenario  in  any 
particular  application.  For  example,  in  the  examples  of  the  following 
section,  we  applied  the  above  computational  procedure  to  Bonder’s 
constant -speed  attack  model  for  which  the  tactical  scenario  is  only 
defined  for  0  <  f  <  =  r0/v.  In  these  examples,  when  the  X  force  is 

not  annihilated  within  this  given  time  t^,  we  calculated  the  final  X 
force  level  by  (6.1)  with  the  help  of  our  tabulations. 

9.  NUMERICAL  EXAMPLES 

In  the  section  we  examine  a  couple  of  numerical  examples  to  show  how 
our  results  lead  to  insights  about  the  dynamics  of  combat  between  two 
homogeneous  forces.  As  in  Taylor  [1974]  and  Taylor  and  Brown  [1976], 
we  consider  Bonder’s  [1965]  model  (3.2)  for  the  constant-speed  attack 
against  a  static  defensive  position.  We  will  focus  on  the  new  results  of 
this  paper  [in  particular,  the  prediction  of  battle  outcome  from  initial 
conditions  without  explicitly  computing  the  force-level  trajectories]  and 
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will  follow  the  computational  procedure  outlined  in  the  previous  section. 
Here  the  tactical  scenario  is  defined  only  for  0  <  t  <  tm,  since  the 
constant-speed  attack  ends  at  =  r0/v.  Hence,  for  t  >  one  must 
not  use  any  results  from  the  mathematical  model  in  which  the  attrition- 
rate  coefficients  (3.5)  have  been  assumed  to  hold  for  all  time.  From  the 
input  data  given  in  Table  II,  we  compute  the  parameter  values  shown  in 
Table  III.  We  observe  from  Table  8B  of  Taylor  and  Brown  [1977b]  and 
Table  III  above,  the  predicted  agreement  between  T(1  —  a)/r(a)  and  the 
limiting  value  of  T„(x )  as  x  — *  +°°  (see  (5.11))  for  a  =  q  —  %  (recall  (4.5)). 
We  now  consider  two  cases:  (I)  r0  =  2000  meters,  and  (II)  r0  =  1250 
meters.  The  interested  reader  can  find  these  examples  worked  out  in 
even  more  detail  in  Taylor  and  Brown  [1977b,  1977c]. 

When  r0  =  2000  meters  (see  Figure  3  of  Taylor  [1974]),  we  have  Ks  = 
0  and  r0  =  0.  The  maximum  time  that  the  battle  can  last  is  t^,  =  14.91 
minutes,  since  at  this  time  the  attackers  reach  their  final  objective  (i.e 
the  defensive  position).  We  now  consider  the  qualitative  behavior  of  the 

TABLE  II 

Input  Data  fob  Numerical  Examples 

It  =  1,  »  =  2. 

a0  =  0.06.X  casualties/minute/  Y  unit. 

/Jo  =  0.6  V  casualties/minute/X  unit. 

r„  =  ra  =  2000  meters. 

0  =  5  miles/hour. 


H  1,  v  =  2  force-level  trajectory  shown  in  Figure  3  of  Taylor  [1974]. 
Theorem  2  tells  us  that  X  can  be  annihilated  if  and  only  if  xo/y0  <  0.420. 
By  (6.3),  the  annihilation  time  of  the  X  force  is  given  by  T,[r(tax)]  = 
3.544  xo/yo-  For  Xo  =  10,  y0  =  30,  we  have  T„(t0x)  =  1.8122  so  that  from 
Table  8A  of  Taylor  and  Brown  [1977b]  (using  linear  interpolation)  we 
obtain  t„x  =  1.009.  Hence,  (4.6)  yields  tBx  =  14.24  minutes  and  rx  =  89.8 
meters.  Further  results  are  given  in  Table  IV. 

When  r„  =  1250  meters  (see  Figure  3  of  Taylor  and  Brown  [1976]),  we 
have  Ks  =  5.5923  meters,  r0  =  0.0975,  and  =  9.32  minutes.  In  this 
case  (again,  for  n  =  1,  v  =  2),  X  can  be  annihilated  if  and  only  if  Xo/yo  < 
0.382  with  (from  (6.3))  the  annihilation  time  of  the  X  force  given  by 
T,(t„x)  =  (3.565uo  +  0.223)/(0.156uo  +  1.004),  where  Uo  =  Xo/y0.  Some 
further  numerical  results  are  given  in  Table  V.  Again,  these  parametric 
results  should  be  contrasted  with  the  single  *i  =  1,  v  *  2  force-level 
trajectory  shown  in  Figure  3  of  Taylor  and  Brown  [1976]. 

10.  FINAL  REMARKS 

In  Sections  6  and  9,  we  have  seen  how  our  new  definition  of  power 
Lanchester  functions  (guided  by  the  general  requirements  for  GLF  given 
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in  Taylor  and  Brown  [1977a])  allows  one  to  conveniently  obtain  much 
valuable  information  about  the  model  (1.1)  with  attrition -rate  coefficients 
(3.5)  without  explicitly  computing  the  entire  force-level  trajectories.  In  his 
well-known  survey  paper  on  the  Lanchester  theory  of  combat,  Dolansky 
suggested  the  development  of  such  outcome-predicting  relations  without 
solving  in  detail  and/or  computing  force-level  trajectories  as  one  of 
several  problems  for  further  research.  Our  Theorem  2  is  a  step  toward 
resolving  this  problem  (see  also  Taylor  and  Parry;  Taylor  and  Comstock; 
Taylor  [1979a]).  Previously,  one  was  limited  to  being  able  to  compute 
only  force-level  trajectories,  but  now  we  can  tell  who  is  going  to  be 
annihilated  (and  when)  without  explicitly  computing  the  trajectories. 

We  have  answered  questions  about  qualitative  model  behavior  (e.g. 
force  annihilation),  not  only  for  specific  values  of,  for  example,  initial 
force  levels,  but  also  for  the  entire  possible  range  of  values  for  the  initial 
force  ratio  (i.e.  parametric  analysis  of  model  behavior).  The  results  of 
this  paper  may  be  used  for  other  parametric  analyses  (see  Bonder  [1971] 
for  a  lucid  discussion  of  the  importance  of  such  analyses),  e.g.  parametric 

TABLE  III 

Parameter  Values  for  Numerical  Examples 

fc.  •  4.0233  X  10-3  X  casualties/! minute)'/!'  unit. 
k»  "  2.6979  x  10-3  Y  casualties/(minute)7X  unit, 
p  =  2/6,  q  =  3/5. 
r(p)/r(9)  -  1.48951. 

Xu  .  0. 


dependence  of  battle  outcome  on  attrition- rate  coefficients.  Thus,  our 
new  results  now  allow  one  to  develop  important  insights  into  the  dynam¬ 
ics  of  combat  between  two  homogeneous  forces  with  temporal  variations 
in  fire  effectiveness.  With  the  availability  (Taylor  and  Brown  [1977b, 
1977c])  of  tabulations  of  the  LCS  functions,  one  can  now  analyze  combat 
modeled  by  the  power  attrition-rate  coefficients  (3.5)  with  somewhat  the 
same  facility  as  he  can  for  the  constant-coefficient  case  and  thus  aid  in 
parameter  analyses. 

In  his  classic  paper,  Lanchester  [1914]  considered  constant  fire  effec¬ 
tiveness  for  individual  firers  and  deduced  his  famous  square  law 

0|*o2  -  x2(t)|  =  a|y02  -  y2(0|,  (10.1) 

where  a  and  0  denote  the  constant  attrition -rate  coefficients.  It  follows 
from  (10.1)  that  (provided  there  is  no  “time  limit”  for  combat) 

X  will  be  annihilated  if  and  only  if  Xo/yo  <  Va//3.  ( 10.2) 

Thus,  we  see  that  equality  of  Lanchester-type  fighting  strengths  depends 
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TABLE  IV 

Annihilation  or  the  X  Force  as  a  Function  or  the  Initial  Force  Ratio 
tor  r„  =  2000  Meters 


(Wyo) 

taX(Minutes) 

r„x(  Meters) 

0.333 

14.24 

89.8 

0.250 

11.61 

443.2 

0.200 

10.19 

633.2 

on  two  parameters:  (I)  initial  force  ratio,  and  (II)  relative  effectiveness. 
When  the  timing  of  military  actions  is  also  considered,  we  add  a  third 
parameter,  the  intensity  of  combat  =  •Saji,  to  this  list  of  significant 
combat  parameters.  No  such  simple  relation  like  the  square  law  (10.1), 
which  yields  (10.2),  holds  in  general  for  variable  attrition-rate  coeffi¬ 
cients.  However,  by  transforming  the  independent  variable  t  to  normalize 
the  battle’s  time  scale  by  the  intensity  of  combat,  we  found  that  the 
course  of  such  variable-coefficient  combat  depends  on  only  the  two 
weapon-system  parameters:  (I)  relative  fire  effectiveness,  R{t)  =  a(t)/ 
b(t),  and  (II)  intensity  of  combat,  I{t)  =  •Ja(t)b(t).  This  way  of  viewing 
the  attrition-rate  coefficients  a(t)  and  b(t)  is  both  intuitively  appealing 
and  also  important  because  under  some  circumstances  relative  fire  effec¬ 
tiveness  (i.e.  only  one  parameter)  plays  the  major  role  in  determining 
battle  outcome  (e.g.  when  a(t)/b(t)  =  constant,  the  intensity  of  combat 
does  not  influence  the  outcome  of  battle  (provided  that  there  is  no  time 
limitation)).  (See  also  Taylor  (1980b].)  Moreover,  we  did  extend  (10.2) 
to  combat  modeled  with  the  power  attrition -rate  coefficients  with  “no 
offset”  (3.5)  (see  Theorem  2).  This  is  the  first  time  that  such  a  generali¬ 
zation  of  the  square  law  has  been  obtained  for  the  variable-coefficient 
Lanchester-type  model  (1.1)  with  a(t)/b(t)  #  constant.  We  observe  that 
for  Ks  >  0,  this  “exact"  outcome-prediction  relation  (i.e.  necessary  and 
sufficient  condition  for  force  annihilation)  involves  higher  transcendental 
functions  (here,  the  LCS  functions)  and  is  complementary  to  the  suffi¬ 
cient  condition  (involving  only  elementary  functions)  given  by  Taylor 
and  Parry  for  Ks  >  0. 

Work  by  Bonder  [1965,  1967,  1970],  Clark  [1969],  Barfoot  [1969],  and 


TABLE  V 

Annihilation  or  the  X  Force  as  a  Function  or  the  Initial  Force  Ratio 
rOR  r*  =  1250  METERS 


(Wyb) 

I„x(Minutes) 

r«x(Meters) 

0.333 

10.63 

_ a 

0.250 

7.56 

235.9 

0.200 

6.17 

422.8 

*  u  =  9.32  minutes  and  xy  = 

x<r  -  0)  =  1.35. 
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Bonder  and  Farrell  on  the  prediction  of  Lanchester  attrition-rate  coef¬ 
ficients  (see  Taylor  and  Brown  [1976]  for  further  discussion  and  refer¬ 
ences)  has  generated  interest  in  variable-coefficient  Lanchester-type 
models.  Interest  in  the  power  attrition-rate  coefficients  with  “no  offset” 
(3.5)  is  provided  by  S.  Bonder’s  [1965]  model  (3.2)  and  his  examination 
of  predicted  attrition  rates  for  various  weapon  systems  (see  pp.  196-200 
of  Bonder  and  Farrell).  However  useful  our  results  may  be  in  their  own 
right,  they  have  far  greater  import:  (I)  they  are  a  model  for  the  treatment 
of  other  Lanchester  functions  and  their  tabulations,  and  (II)  they  may 
be  used  in  the  numerical  determination  of  the  parity-condition  parameter 
Q *  for  related  attrition -rate  coefficients  (e.g.  (3.1)  with  K0  >  0).  In  Taylor 
and  Brown  [1978],  we  show  how  our  tabulations  of  the  LCS  functions 
play  a  key  role  in  the  numerical  determination  of  the  parity  condition 
parameter  Q*  for  the  general  power  attrition -rate  coefficients  (3.1)  with 
positive  “offset”  (i.e.  Ko  >  0). 

We  have  extended  our  mathematical  theory  (Taylor  and  Brown 
[1976])  of  variable-coefficient  Lanchester-type  equations  of  “modern 
warfare”  for  combat  between  two  homogeneous  forces  in  order  to  be  able 
to  more  thoroughly  analyze  such  models  (see  also  Taylor  and  Brown 
[1977a]).  The  classic  ordinary-differential-equation  theories  (e.g.  see 
Hille  [1969])  were  inadequate  to  supply  all  the  answers  sought  about 
such  combat  models  (Taylor  [1979b]).  The  mathematical  theory  of  the 
model  (1.1)  with  coefficients  (3.5)  is  now  nearly  as  complete  as  that  of 
the  constant-coefficient  model.  Such  results  as  we  have  given  here  are 
very  useful  for  understanding  the  dynamics  of  combat,  i.e.  how  the  trading 
of  casualties  will  be  projected  over  time.  H.  K.  Weiss  [1959]  has  empha¬ 
sized  that  such  a  simplified  combat  model  is  particularly  valuable  when 
it  leads  to  a  clearer  understanding  of  such  significant  relationships  that 
would  tend  to  be  obscured  in  a  more  complex  model. 
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